Block Diagonal and Schur Complement Preconditioners for Block-Toeplitz Systems with Small Size Blocks
نویسندگان
چکیده
Abstract. In this paper we consider the solution of Hermitian positive definite block-Toeplitz systems with small size blocks. We propose and study block diagonal and Schur complement preconditioners for such block-Toeplitz matrices. We show that for some block-Toeplitz matrices, the spectra of the preconditioned matrices are uniformly bounded except for a fixed number of outliers and this fixed number depends only on the size of block. Hence conjugate gradient type methods, when applied to solving these preconditioned block-Toeplitz systems with small size blocks, converge very fast. Recursive computation of such block diagonal and Schur complement preconditioners are considered by using the nice matrix representation of the inverse of a block-Toeplitz matrix. Applications to block-Toeplitz systems arising from least squares filtering problems and queueing networks are presented. Numerical examples are given to demonstrate the effectiveness of the proposed method.
منابع مشابه
A new approximation of the Schur complement in preconditioners for PDE-constrained optimization
Saddle point systems arise widely in optimization problems with constraints. The utility of Schur complement approximation is now broadly appreciated in the context of solving such saddle point systems by iteration. In this short manuscript, we present a new Schur complement approximation for PDE constrained optimization, an important class of these problems. Block diagonal and block triangular...
متن کاملBlock Diagonal Preconditioners for the Schur Complement Method Block Diagonal Preconditioners for the Schur Complement Method
We present numerical methods for solving systems of linear equations originated from the discretisation of two-dimensional elliptic partial diierential equations. We are interested in diierential equations that describe heterogeneous and anisotropic phenomena. We use a nonoverlapping domain decomposition method for solving the linear systems. We describe new local preconditioners for the interf...
متن کاملA Direct Elliptic Solver Based on Hierarchically Low-rank Schur Complements
Cyclic reduction was conceived for the solution of tridiagonal linear systems, such as the one-dimensional Poisson equation [11]. Generalized to higher dimensions, it is known as block cyclic reduction (BCR) [4]. It can be used for general (block) Toeplitz and (block) tridiagonal linear systems; however, it is not competitive for large problems, because its arithmetic complexity grows superline...
متن کاملDevelopments in Preconditioned Iterative Methods with Application to Glacial Isostatic Adjustment Mo- dels
This study examines the block lower-triangular preconditioner with element-wise Schur complement as the lower diagonal block applied on matrices arising from an application in geophysics. The element-wise Schur complement is a special approximation of the exact Schur complement that can be constructed in the finite element framework. The preconditioner, the exact Schur complement and the elemen...
متن کاملRobust Preconditioned Iterative So- lution Methods for Large-scale Non- symmetric Problems
We study robust, preconditioned, iterative solution methods for largescale linear systems of equations, arising from different applications in geophysics and geotechnics. The first type of linear systems studied here, which are dense, arise from a boundary element type of discretization of crack propagation in brittle material. Numerical experiment show that simple algebraic preconditioning str...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- SIAM J. Matrix Analysis Applications
دوره 29 شماره
صفحات -
تاریخ انتشار 2007